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Spectroscopy has played the key role in revealing, and thereby understanding, the structure of 
atoms and molecules. A central drive in this field is the pursuit of higher precision and accuracy so 
that ever more subtle effects might be discovered. Here, we report on laser absorption spectroscopy 
that operates at the conventional quantum limit imposed by photon shot-noise. Furthermore, we 
achieve this limit without compromising the accuracy of the measurement. We demonstrate these 
properties by recording an absorption profile of cesium vapor at the 2 parts-per-million level. The 
extremely high signal-to-noise ratio allows us to directly observe the homogeneous lineshape compo- 
nent of the spectral profile, even while in the presence of Doppler broadening that is a factor of 100 
times wider. We can do this because we can precisely measure the spectral profile at a frequency 
detuning more than 200 natural line widths from the line center. We use the power of this tool to 
demonstrate direct measurements of a low-intensity optically-induced broadening process that is 
quite distinct from the well-known power broadening phenomenon. 



Researchers performing absorption spectroscopy have 
continually sought to improve the precision with which 
atomic and molecular lines can be measured. As instru- 
ment precision improved, comparisons of the observed 
spectra to contemporary theory often led to the discovery 
of subtle structures that were previously unresolvable. 
Exemplary cases include the observation, and consequen- 
tial understanding, of the Stark[l and Lamb shifts [2 . 
Attaining accuracy as well as precision in modern fre- 
quency metrology now rests on the detailed removal of 
these and many similar systematic frequency effects at 
the level of a few parts in 10^^[3l[4]. Analogously, recent 
attainment of astonishing frequency precision and purity 
in the modes of an optical frequency comb are offering 
new possibilities for laboratory-scale searches for varia- 
tions in the values of fundamental constants[5^ Here, 
it is again demonstrated that achieving incredible pre- 
cision is only half the story; accuracy at the same level 
is also required, which demands stringent understand- 
ing of all systematic factors in the model for the ob- 
served spectroscopic lineshape [6 . Progress in precision 
laser absorption spectroscopy has thus become an en- 
abling technology of key importance in a wide variety 
of fields, which include frequency metrology [71 [8], di- 
rect optical frequency comb spectroscopy p) [TO] , primary 
thermometry [TTHIl], trace gas detection J5l [16], preci- 



sion lineshape measurements [TTl [TSj and precision line- 
shape analysis [19 . 

The conventional limit to precision in laser absorption 
spectroscopy (LAS) is set by shot-noise in the probing 
light, though, in practice, it is difficult to achieve this 
limit as technical noise and instrumental limitations usu- 
ally intervene. These technical limitations are partic- 
ularly acute for LAS because it is inherently a bright 
field[20^ measurement i.e. the signal of interest is a 
difference between two large signals: one proportional 
to the incident power, and the second proportional to 
the transmitted power. If one wishes to build a high- 
sensitivity and high-accuracy tool for probing optical ab- 
sorbance then there are two key challenges: first, small 
signals are easily obscured by technical noise sources such 
as laser amplitude instability, which are typically much 
larger than fundamental sources; and second, the resolu- 
tion of the measurement (the small differential signal) is 
compromised because most of the dynamic range of the 
sensor is consumed by the need to measure the large av- 
erage value of the signal. Ingenious differential-detection 
designs, such as the noise eater described in Ref. [21], 
can be used to reach the shot-noise limit in the presence 
of amplitude noise that can up to 70 dB larger through a 
high degree of common mode rejection. This type of bal- 
anced detector was recently demonstrated in conjunction 



with an optical frequency comb to give high-precision 
broadband measurements at the shot-noise hmit[22 . 

Although the precision of such techniques is extremely 
impressive, it is important to note that this often comes 
at the expense of accuracy because information about the 
average power (common mode signal) is either not mea- 
sured, or because the design of the differential detector 
does not adequately ensure linearity. The latter prob- 
lem is typically the case for auto-balancing photocurrent 
attenuators [21] [23| [24] . This loss in accuracy is not crit- 
ical if the ultimate goal is to achieve extraordinary sen- 
sitivity; however, in many areas of current interest, ac- 
curacy is just as important as precision [25] |26]. In such 
circumstances the strict adherence of the measured line- 
shape to the assumed model profile is critical to avoid- 
ing systematic errors when line parameters are extracted 
from regression. To this end, other work has previously 
sought, and found, significant non-linearity in the absorp- 
tion depth of alkali metal vapors even in the weak-probe 
regime (/ < L^^ /lQ)[27r.29j. However, these earlier stud- 
ies were unable to resolve the subtle changes in homo- 
geneous linewidth that are at the origin of the lineshape 
changes. 

One brute-force approach to overcome this challenge 
would be to use traditional LAS but digitize the two sig- 
nals of interest with an extraordinarily large dynamic 
range device (a possibility canvased in Ref. [22]). How- 
ever, as we demonstrate in this paper, it is possible with 
a simple rearrangement of existing equipment, to obtain 
a linear differential LAS technique that allows measure- 
ments of optical power ratios at the few parts-per- million 
level with precision limited only by shot-noise on the 
light. 

We make use of this high-resolution spectrometer to 
probe the Dl transition of cesium (Cs) and report an 
unexpected source of optically-induced modification to 
the homogeneous lineshape component in the absorption 
profile. Whilst it is well established that the homoge- 
neous linewidth can be substantially modified at inten- 
sities higher than the saturation intensity ^Ql |3Tj, ^sat, 
we directly observe that the linewidth has doubled at in- 
tensities as low as /gat/ 10- This effect is an unforeseen 
consequence of optical pumping in the alkali vapours. 

Our high-accuracy, high-precision technique can find 
application in lineshape analysis to extract collisionally- 
induced perturbation parameters [ITOl l32]. verification of 
time-dependent optical pumping models [33]. and studies 
of plasma condition by lineshape analvsis[34][35] 



RESULTS 

In conventional LAS detection, the incident and trans- 
mitted light is independently detected (see inset on Fig. 
Tl) to deliver a photocurrent of / and T respectively [13j 
31j . Any common- mode laser intensity fluctuations are 
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FIG. 1: (Color Online) Quasi-differential detection scheme for 
highly linear, large dynamic range bright field measurements 
of spectral lines. A frequency-stabilized laser is split at the 
Wollaston prism. The reference beam is detected on photo- 
diode /. The probe beam passes through the Cs cell inside a 
thermally and magnetically controlled environment, and the 
transmitted light is detected on photodiode T. Inset: Tradi- 
tional LAS detection technique with sub-optimal resolution 
limited by instrumental dynamic range. The optical setup is 
unchanged. 



removed by constructing the ratio T/I. This procedure 
intrinsically compromises the resolution of an absorption 
measurement as most of the dynamic range of the pho- 
tocurrent measurement device is used to deliver solely in- 
formation on the average photocurrent. In contrast, our 
technique is based on effectively extending the dynamic 
range of the photocurrent measurement by making two 
different measurements: first, a high-gain measurement 
of the difference signal I — T which is centred on zero, 
and a separate measurement of the average component 
(/). The ratio T/I can then be reconstructed in software 
using the identity 



= 1 



I-T 

I ' 



(1) 



The increase in dynamic range of this new technique 
over the traditional method is equal to the inverse of 
the depth of the absorption signal (i.e. in our case 
^ 1/0.05 = 20). This condition is set by the maximum 
increase in gain allowable on the detection of / — T over 
that of / by itself. It is this careful differential approach 
that has given us a linear, high dynamic-range ratiomet- 
ric measurement of the optical powers. We have verified 
that non-linearity of the detection chain contributes on 
the order of 100 parts-per-million detection gain varia- 
tion between 95 and 100 % of the typical operating opti- 
cal powers, (a few microwatts) which we can correct to 
first order to the level of the uncertainty of our linear- 
ity verification (a few parts-per-million) (see Methods). 
We expect that the photocurrent ratio (T/I) is an ac- 
curate measurement of the transmission P/Pq to within 
this limit. 

Our measurements are made at probe intensities be- 
tween 10~^/sat to 10~^/sat of the Saturation intensity. 



where /gat = 2.5mW/cm^j36]. At such low probe inten- 
sities, the transmitted probe power P{f) can be modeled 
using the Beer-Lambert law[30l[3l], 



P(/) = Po(/)exp[-<T], 



(2) 



where Po{f) is the probe power incident on the atoms at 
each optical frequency / and a is the (frequency depen- 
dent) optical depth. In this same low intensity limit, we 
can write a = aV{f — /o), where a is the on-resonance 
optical depth; /o is the optical frequency of the transi- 
tion and V{f) is a Voigt function [171 E] that describes 
the real-part of the complex susceptibility of the opti- 
cal resonance for a thermal ensemble of atoms. The 
Voigt function is a convolution of a Lorentzian com- 
ponent of the lineshape (homogeneous linewidth) and 
a Gaussian component (inhomogenous Doppler broad- 
ening). The Lorentzian component can be defined as 
Hf) = [1 + (//r)^]"^) where F is the half- width at 
half maximum (HWHM) bandwidth, while the Gaus- 
sian component, G{f) = exp(— //(^^y) has a width of 
5uj = \/ {2kT) / {mc^) where k is the Boltzmann constant, 
T is the temperature of the vapor, m is the atomic mass 
and c is the speed of light. Broadening effects from tran- 
sit time and the probing laser are negligible in this ex- 
periment. Power broadening will only amount to an ap- 
parent 5% change in the homogenous linewidth over the 
probe intensity range [30 . We thus expect the Lorentzian 
component to be nearly constant with a width set by 
he natural lifetime of the transition, i.e. F = 2.29 MHz. 
All of our measurements were made at T = 279 K so 
that all of the recorded spectra should also have a fixed 
Gaussian component with the same characteristic width 
(-207 MHz). 

A typical transmission spectrum (T/I) for which / = 
0.15/sat and T = 279.29 ± 0.015 K is shown in the lower 
panel of Fig. [2] A least-squares fit was performed to 
the spectrum and the residual from a fit to the Eqn. 
[2] model is shown in the upper panel of Fig. [2] The 
residuals are featureless and have a white spectrum with 
an root-mean-square (RMS) spread of 2ppm [44]. This is 
entirely consistent with the expected photon shot-noise 
at this probe power (6 /iW). 

To demonstrate explicitly that the measurement noise 
is associated with shot-noise we monitored the fiuctua- 
tions in the measurement as the probe laser was held at 
three different frequencies corresponding to a) completely 
out of resonance with the Cs, b) at a point halfway down 
the absorption curve, and c) tuned to the centre of reso- 
nance. Each of these measurements should yield similar 
results if frequency noise of the laser source is negligible 
(since the transmitted power is similar in all cases). Fig. 
[3] (top) shows the measured fiuctuations in each of these 
signals over various time scales. We display this in terms 
of the Allan Variance which is a standard and stable mea- 
sure of fiuctuations in an experimentally obtained time 
series 138^. The measurement noise is seen to be nearly 
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FIG. 2: (Color Online) Lower panel: A typical spectrum of 
a Dl Cs absorption line. Open circles are the data, and the 
solid blue is the fit. Upper panel: Residuals to the fit func- 
tion in Eqn. [2] and is consistent with a shot-noise limited 
measurement. 



identical in all cases confirming that the laser frequency 
noise is negligible at the current measurement precision. 
The slope of each result between 40 ms and 0.2 s of in- 
tegration time indicates that the fiuctuation spectrum is 
white. For longer timescales, the measurement is dom- 
inated by fiicker-noise arising in the electronic compo- 
nents of the measurement system and possible drifts in 
the alignment of stray etalons[39 . Fig. [sFb) shows the 
measurement noise at 40 ms integration time as a func- 
tion of probe power. Here, we see a distinct change in 
slope around 3.5/iW incident power from 1/P to 1/VP 
behaviour. We also display the calculated Johnson noise 
of the photodetector readout resistor (blue) and the cal- 
culated shot-noise of the light from first principles (green 
curve). It can be seen that the dependence and magni- 
tude of the noise at low powers is consistent with Johnson 
noise, while the shape and magnitude is consistent with 
photon shot-noise at optical powers higher than 3.5 /iW. 
The large dynamic range and high SNR of our ap- 
proach means that we can directly reveal subtle and im- 
portant changes in the absorption spectra, which have 
not been observed previously. Fig.|4]shows the frequency- 
dependent optical depth, a, extracted from the data us- 
ing Eqn. [2J We display this depth for a set of probe 
powers ranging from ~ 10~^ — 10~^/sat- One observes 
a strong on-resonance Gaussian profile associated with 
Doppler broadening of the vapour. In addition, the high 
signal to noise ratio in our experiment enables us to ob- 
serve the Lorentzian wings that make up a "pedestal" 
of the Voigt profile at large detunings > lOOF. These 
are easily distinguished: the Gaussian component decays 
rapidly at high frequency detunings until the Lorentzian 
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FIG. 3: Noise in the measurement of optical transmission 
relative to the off-resonance value of near unity. (Top) The 
relative noise is white on timescales between 40 ms to 200 ms, 
but becomes dominated by flicker noise at longer times. There 
is almost no additional noise induced by frequency instability 
of the probe laser. (Bottom) The relative noise in the mea- 
surement of the optical transmission ratio at an integration 
time of 40 ms. The measurements become shot-noise limited 
at probe powers greater than 3 /iW. 



wings become dominant [40| . 

Fig. [4] shows a surprising, and at first sight quite un- 
expected, feature: it can be clearly seen that magnitude 
of the Lorentzian component is dependent on the probe 
intensity. This is even more evident on the inset on this 
figure where we re-plot the data in terms of its appar- 
ent Lorentzian width (F): this increases from 2 MHz to 
6 MHz as the probe intensity increases from 10~^/sat to 
10~^/sat whereas we should have expected that this was 
constant within 5% over this intensity range [36 . We 
have recently shown that this type of behaviour can arise 
out of a complex interplay between optical pumping and 
atom dynamics as an atom traverses the probe beam [37 . 
This interaction modifies the underlying Lorentzian line- 
shape to give an overall lineshape that is still close to 
a Voigt function but with a modified Lorentzian width 
parameter. 
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FIG. 4: (Color Online) Demonstration of the breakdown in 
the Voigt profile at low intensities, (top) Optical depth versus 
detuning for various probe intensities. The inset shows the 
fitted Lorentzian width as a function of probe power; the fact 
that this is not constant over the range of probe intensities 
is a signature of the breakdown of the Voigt spectral line 
shape [37] . The reference profile which is expected in the limit 
of zero probe intensity consists of a Gaussian component of 
207.902 MHz and a Lorentzian component of 2.2875 MHz. 



DISCUSSION 



We have demonstrated a new approach to laser absorp- 
tion spectroscopy that circumvents the dynamic range 
limitations of the standard technique. This is achieved by 
using a high-precision differential measurement scheme 
that carefully preserves the measurement linearity while 
simultaneously obtaining scale information in order to 
deliver accuracy. This has permitted the reconstruc- 
tion of Cs absorption spectra with a SNR limited by 
shot-noise on the probe source (2ppm at 6/iW in a 
40ms bandwidth). By varying the probe power between 
1/100 to 1/10 of the saturation intensity, we have ob- 
served an apparent trebling of the underlying homoge- 
neous linewidth of the Cs Dl transition. This observation 
shows the power of the new technique as we are making 
precision measurements of the natural linewidth compo- 
nent of the voigt profile in the wings of the spectral line at 
frequency detuning of more than 200 natural linewdiths 
away from the line center. This work heralds a potential 
leap in ultra-high precision absorption spectroscopy of 
gases where we can now overcome previously overlooked 
physical phenomena that arise from time-dependent op- 
tical pumping. We anticipate these techniques will be 
critical in future high precision spectroscopy research. 




FIG. 5: (Color Online) Optical setup. The inset showns a 
magnified view of the output window of the Cs cell. There is 
a small volume overlap, shaded orange, between the incoming 
probe beam and its weak retroreflection from the angled, but 
uncoated cell window. ^ = 15°. 



METHODS 
Optical setup 

A detailed schematic of the experimental setup is 
shown in Fig. |5j One extended cavity diode laser 
(ECDL), labeled reference^ was frequency-stabilized to 
a saturated-absorption hyperfine spectroscopic feature of 
Cs at 894.578 nm in the Dl manifold, using a standard 
Pound-Drever-Hall (PDH) laser- locking technique [41 . A 
second ECDL, labeled probe, was used to interrogate the 
cesium transitions at 894.580 nm, which are 1.17 GHz de- 
tuned from the reference transition. This laser was fre- 
quency locked using the PDH technique to the funda- 
mental transverse mode of an optical cavity (OC) (1 MHz 
bandwidth; 18 GHz free spectral range). A beat-note be- 
tween the two lasers was formed on a high speed pho- 
todiode (DB), whose output was frequency divided by a 
factor of n = 20480, and locked to a tunable radio fre- 
quency oscillator {/rf) using a frequency- locked loop to 
stabilize the frequency probe at the 2 kHz level over 1 
to 10 s. By manipulating fjip the probe laser frequency 
could be tuned to an arbitrary frequency offset from the 
reference laser. A frequency counter was used to mon- 
itor the (divided) beat-note frequency, fbeat^ which was 
used to accurately reconstruct the frequency axis for the 
transition of interest. 

Light transmitted by the optical cavity, which was 
mostly free of spontaneous emission from the ECDL, was 
passed through an acousto-optic modulator (AOM). This 
AOM performed two tasks: an optical chopper (90% 
amplitude modulation at fmod = 1.57 kHz), and as a 



variable optical attentuator to deliver high-bandwidth 
power stabilization of the diffracted beam. Th diffracted 
beam is shifted in frequency by the AOM's drive fre- 
quency of 175 MHz, f drive- The zero-order undiffr acted 
beam was used to compare the probe laser frequency to 
a stabilised optical frequency comb (not shown). The 
diffracted beam was coupled into a length of single-mode 
optical fiber (OF) that eliminated pointing fluctuations 
in the light. The resulting spatially and spectroscopi- 
cally pure probe beam was then polarized and split using 
Clan- Taylor (GT) and Wollaston (WP) prisms respec- 
tively. Light from one arm was detected immediately, 
while light from the other arm was directed through the 
cesium cell before detection (on DI and DT respectively) . 

The detectors were reverse-biased silicon photodiodes, 
with matched 47 kl^ load resistors. The voltage across 
each load resistor was synchronously detected using a 
lock-in amplifier (referenced to fmod)^ and recorded on a 
digital multimeter (DMM). To circumvent the dynamic 
range limit imposed by the lock-in amplifier (Stanford 
Research Systems SR830), which was 10 times worse than 
the shot-noise limit, a low- noise analog subtraction stage 
was used to create a differential signal that was close to 
zero when the laser was off-resonance. Lock-in ampli- 
fier 2 (LI 2) was set to demodulate the differential signal 
Vsig — Vref- The advantage of this technique is its abil- 
ity to measure small transmission deviations from unity 
without loss of dynamic range. In addition, the incident 
photodiode (DI) output was detected by a third lock-in 
amplifier to produce a signal for optical power stabiliza- 
tion. The output of this third lockin was compared to 
a reference voltage and then fed back to the AOM via 
a voltage controlled attenuator (VGA) to stabilize the 
optical power of the probe beam. 



The cesium cell under test was encased in a 15 kg cylin- 
drical copper block. Its temperature was monitored using 
a using two platinum resistance thermometers (PRTs) 
calibrated by their manufacturer with an uncertainty of 
30 mK. One PRT was placed at the centre of the block 
(shown in Fig. Isl), 5 mm from the mid-point of the cell; 
the other was located at the extreme end of the block 
(not shown). A temperature-controlled chiller sent re- 
frigerated liquid through a copper pipe coiled around the 
block to maintain a temperature of approximately 276 K, 
with stability of 3 mK over one minute and long term sta- 
bility of 15 mK. At this temperature, the cesium vapor in 
the cell was optically thin. The maximum gradient ob- 
served between the PRTs was 30 mK, which was at the 
same level as their combined calibration uncertainty. The 
copper block and cooling coils were enclosed within a mu- 
metal magnetic shield to minimize Zeeman broadening 
and mounted on sorbothane pads for passive vibration 
isolation. 



Detection Linearity 

The linearity of the entire detection scheme (reverse- 
biased photodiode and load resistor; lock-in amplifier; 
and digital multimeter) was determined using a mod- 
ified, synchronous version of the method discussed by 
Shin [43 . Two lasers of equal optical power were com- 
bined and passed through an AOM (to produce AM mod- 
ulation identical to that discussed above) before falling 
on the photodiode. Each beam could be independently 
switched on or off using a mechanical chopper. A third 
laser, with optical power 20 times that of each chopped 
laser, was also passed through the AOM and onto the 
detector to introduce a signal replicating our operating 
point. The two low-power beams were continuously cy- 
cled through the four possible states (on:off, on:on, off:on, 
offioff) and the resulting synchronously detected output 
voltages were recorded. The figure of merit for this lin- 
earity measurement is the relative error, given by: 




Number of Spectra 

FIG. 6: (Color Online) Relative Allan deviations for selected 
fit parameters. The solid lines are extracted from measured 
spectra, whilst the dashed traces are from simulations using 
noise characteristics consistent with that shown in Fig. [3k. 



{Vi + V2)- V12 



(3) 



where Vi, V2 and F12 are the recorded output voltages 
for the first beam only, the second beam only and both 
beams together, respectively. The relative error can be 
calculated for each complete cycle with the fourth state 
(both beams off) used to remove unwanted offsets. Av- 
eraging the results of many such measurements reduces 
uncertainty in the relative error to arbitrarily small levels. 
By measuring the relative error as a function of optical 
power (or more practically, output voltage) the resulting 
first-order non-linearity of the detection scheme can be 
deduced and corrected. 



Data reduction 

In the least-squares fits of the measured transmission 
data, a small quadratic power dependence arising from 
broad and weak etalons formed by the surfaces of each 
uncoated cell window. This is modeled by allowing the 
input power in Eqn. 0to take the form Po(/) = ^/^ + 
Bf + C . The parameters A, B and C are adjustable 
fit parameters. To display the underlying Voigt profiles 
(FigQ, the parameters A^B^C and a are substituted 
with their fitted values into 



n/-/o) = -HP(/)/Po(/)]/a. 



Reproducibility of fit parameters 



(4) 



We recorded a set of 60 consecutive scans of the Cs 
spectrum at a probe power of 6/iW where the measure- 
ment is shot-noise limited. Fig. [6] displays the relative 



Allan deviation of a selection of fit parameters in Eqn. [4] 
(C the off-resonance ratio, a the absorption depth, F the 
Lorentzian half- width and the Gaussian half- width). The 
dashed traces are the expected parameter noises given 
the noise characteristics of the transmission noise (Fig. 
|3^) and demonstrates that stability of the fitted parame- 
ters is consistent with the shot-noise limit in combina- 
tion with the fiicker- noise at longer time scales. The 
Lorentzian width estimation is 100-fold noisier than the 
Gaussian because they are highly fit correlated parame- 
ters; a unit change in the absolute width of one compo- 
nent induces a near unit change in the other. Therefore, 
the relative fiuctuations are in the ratio of the mean val- 
ues of these parameters (~ 100 : 1). Similarly, the pa- 
rameters associated with amplitudes, the optical depth a 
and the optical split ratio C, have similar absolute devi- 
ations from their mean values. However, the mean value 
of the depth is ~ 1/20 of the background level due to 
the chosen cell temperature. The fiicker-noise prevent 
the uncertainty in the fitted parameters decreases with 
ensemble averaging after a few scans. 
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